A simple model for predicting soil temperature in snow-covered and seasonally frozen soil: model description and testing

Microbial processes in soil are moisture, nutrient and temperature dependent and, consequently, accurate calculation of soil temperature is important for modelling nitrogen processes. Microbial activity in soil occurs even at sub-zero temperatures so that, in northern latitudes, a method to calculate soil temperature under snow cover and in frozen soils is required. This paper describes a new and simple model to calculate daily values for soil temperature at various depths in both frozen and unfrozen soils. The model requires four parameters average soil thermal conductivity, specific beat capacity of soil, specific heat capacity due to freezing and thawing and an empirical snow parameter. Precipitation, air temperature and snow depth (measured or calculated) are needed as input variables. The proposed model was applied to five sites in different parts of Finland representing different climates and soil types. Observed soil temperatures at depths of 20 and 50 cm (September 1981-August 1990) were used for model calibration. The calibrated model was then tested using observed soil temperatures from September 1990 to August 2001. R-2-values of the calibration period varied between 0.87 and 0.96 at a depth of 20 cm and between 0.78 and 0.97 at 50 cm. R-2 -values of the testing period were between 0.87 and 0.94 at a depth of 20cm. and between 0.80 and 0.98 at 50cm. Thus, despite the simplifications made, the model was able to simulate soil temperature at these study sites. This simple model simulates soil temperature well in the uppermost soil layers where most of the nitrogen processes occur. The small number of parameters required means, that the model is suitable for addition to catchment scale models.

The diurnal course of soil moisture as measured by various dielectric sensors: Effects of soil temperature and the implications for evaporation estimates

Soil moisture content, theta, of a bare and vegetated UK gravelly sandy loam soil (in situ and repacked in small lysimeters) was measured using various dielectric instruments (single-sensor ThetaProbes, multi-sensor Profile Probes, and Aquaflex Sensors), at depths ranging between 0.03 and I m, during the summers of 2001 (in situ soil) and 2002 (mini-lysimeters). Half-hourly values of evaporation, E, were calculated from diurnal changes in total soil profile water content, using the soil water balance equation. For the bare soil field, Profile Probes and ML2x ThetaProbes indicated a diurnal course of theta that did not concur with typical soil physical observations: surface layer soil moisture content increased from early morning until about midday, after which theta declined, generally until the early evening. The unexpected course of theta was positively correlated to soil temperature, T-s, also at deeper depths. Aquaflex and ML1 ThetaProbe (older models) outputs, however, reflected common observations: 0 increased slightly during the night (capillary rise) and decreased from the morning until late afternoon (as a result of evaporation). For the vegetated plot, the spurious diurnal theta fluctuations were less obvious, because canopy shading resulted in lower amplitudes of T-s. The unrealistic theta profiles measured for the bare and vegetated field sites caused diurnal estimates of E to attain downward daytime and upward night-time values. In the mini-lysimeters, at medium to high moisture contents, theta values measured by (ML2x) ThetaProbes followed a relatively realistic course, and predictions of E from diurnal changes in vertically integrated theta generally compared well with lysimeter estimates of E. However, time courses of theta and E became comparable to those observed for the field plots when the soil in the lysimeters reached relatively low values of theta. Attempts to correct measured theta for fluctuations in T, revealed that no generally applicable formula could be derived. (c) 2005 Elsevier B.V. All rights reserved.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - I. A universal approach to calculate thermal inertia

In this paper we pledge that physically based equations should be combined with remote sensing techniques to enable a more theoretically rigorous estimation of area-average soil heat flux, G. A standard physical equation (i.e. the analytical or exact method) for the estimation of G, in combination with a simple, but theoretically derived, equation for soil thermal inertia (F), provides the basis for a more transparent and readily interpretable method for the estimation of G; without the requirement for in situ instrumentation. Moreover, such an approach ensures a more universally applicable method than those derived from purely empirical studies (employing vegetation indices and albedo, for example). Hence, a new equation for the estimation of Gamma(for homogeneous soils) is discussed in this paper which only requires knowledge of soil type, which is readily obtainable from extant soil databases and surveys, in combination with a coarse estimate of moisture status. This approach can be used to obtain area-averaged estimates of Gamma(and thus G, as explained in paper II) which is important for large-scale energy balance studies that employ aircraft or satellite data. Furthermore, this method also relaxes the instrumental demand for studies at the plot and field scale (no requirement for in situ soil temperature sensors, soil heat flux plates and/or thermal conductivity sensors). In addition, this equation can be incorporated in soil-vegetation-atmosphere-transfer models that use the force restore method to update surface temperatures (such as the well-known ISBA model), to replace the thermal inertia coefficient.

A sensitivity analysis of soil moisture retrieval from the t-w microwave emission model

The potential of the τ-ω model for retrieving the volumetric moisture content of bare and vegetated soil from dual polarisation passive microwave data acquired at single and multiple angles is tested. Measurement error and several additional sources of uncertainty will affect the theoretical retrieval accuracy. These include uncertainty in the soil temperature, the vegetation structure and consequently its microwave singlescattering albedo, and uncertainty in soil microwave emissivity based on its roughness. To test the effects of these uncertainties for simple homogeneous scenes, we attempt to retrieve soil moisture from a number of simulated microwave brightness temperature datasets generated using the τ-ω model. The uncertainties for each influence are estimated and applied to curves generated for typical scenarios, and an inverse model used to retrieve the soil moisture content, vegetation optical depth and soil temperature. The effect of each influence on the theoretical soil moisture retrieval limit is explored, the likelihood of each sensor configuration meeting user requirements is assessed, and the most effective means of improving moisture retrieval indicated.

Effects of soil conditions and drought on egg hatching and larval survival of the clover root weevil (Sitona lepidus)

Soil-dwelling insect herbivores are significant pests in many managed ecosystems. Because eggs and larvae are difficult to observe, mathematical models have been developed to predict life-cycle events occurring in the soil. To date, these models have incorporated very little empirical information about how soil and drought conditions interact to shape these processes. This study investigated how soil temperature (10, 15, 20 and 25 °C), water content (0.02 (air dried), 0.10 and 0.25 g g−1) and pH (5, 7 and 9) interactively affected egg hatching and early larval lifespan of the clover root weevil (Sitona lepidus Gyllenhal, Coleoptera: Curculionidae). Eggs developed over 3.5 times faster at 25 °C compared with 10 °C (hatching after 40.1 and 11.5 days, respectively). The effect of drought on S. lepidus eggs was investigated by exposing eggs to drought conditions before wetting the soil (2–12 days later) at four temperatures. No eggs hatched in dry soil, suggesting that S. lepidus eggs require water to remain viable. Eggs hatched significantly sooner in slightly acidic soil (pH 5) compared with soils with higher pH values. There was also a significant interaction between soil temperature, pH and soil water content. Egg viability was significantly reduced by exposure to drought. When exposed to 2–6 days of drought, egg viability was 80–100% at all temperatures but fell to 50% after 12 days exposure at 10 °C and did not hatch at all at 20 °C and above. Drought exposure also increased hatching time of viable eggs. The effects of soil conditions on unfed larvae were less influential, except for soil temperature which significantly reduced larval longevity by 57% when reared at 25 °C compared with 10 °C (4.1 and 9.7 days, respectively). The effects of soil conditions on S. lepidus eggs and larvae are discussed in the context of global climate change and how such empirically based information could be useful for refining existing mathematical models of these processes.

Using Distributed Temperature Sensing to monitor field scale dynamics of ground surface temperature and related substrate heat flux

We present one of the first studies of the use of Distributed Temperature Sensing (DTS) along fibre-optic cables to purposely monitor spatial and temporal variations in ground surface temperature (GST) and soil temperature, and provide an estimate of the heat flux at the base of the canopy layer and in the soil. Our field site was at a groundwater-fed wet meadow in the Netherlands covered by a canopy layer (between 0-0.5 m thickness) consisting of grass and sedges. At this site, we ran a single cable across the surface in parallel 40 m sections spaced by 2 m, to create a 40×40 m monitoring field for GST. We also buried a short length (≈10 m) of cable to depth of 0.1±0.02 m to measure soil temperature. We monitored the temperature along the entire cable continuously over a two-day period and captured the diurnal course of GST, and how it was affected by rainfall and canopy structure. The diurnal GST range, as observed by the DTS system, varied between 20.94 and 35.08◦C; precipitation events acted to suppress the range of GST. The spatial distribution of GST correlated with canopy vegetation height during both day and night. Using estimates of thermal inertia, combined with a harmonic analysis of GST and soil temperature, substrate and soil-heat fluxes were determined. Our observations demonstrate how the use of DTS shows great promise in better characterising area-average substrate/soil heat flux, their spatiotemporal variability, and how this variability is affected by canopy structure. The DTS system is able to provide a much richer data set than could be obtained from point temperature sensors. Furthermore, substrate heat fluxes derived from GST measurements may be able to provide improved closure of the land surface energy balance in micrometeorological field studies. This will enhance our understanding of how hydrometeorological processes interact with near-surface heat fluxes.

Adaptation of the Integrated Nitrogen Model for Catchments (INCA) to seasonally snow-covered catchments

Testing of the Integrated Nitrogen model for Catchments (INCA) in a wide range of ecosystem types across Europe has shown that the model underestimates N transformation processes to a large extent in northern catchments of Finland and Norway in winter and spring. It is found, and generally assumed, that microbial activity in soils proceeds at low rates at northern latitudes during winter, even at sub-zero temperatures. The INCA model was modified to improve the simulation of N transformation rates in northern catchments, characterised by cold climates and extensive snow accumulation and insulation in winter, by introducing an empirical function to simulate soil temperatures below the seasonal snow pack, and a degree-day model to calculate the depth of the snow pack. The proposed snow-correction factor improved the simulation of soil temperatures at Finnish and Norwegian field sites in winter, although soil temperature was still underestimated during periods with a thin snow cover. Finally, a comparison between the modified INCA version (v. 1.7) and the former version (v. 1.6) was made at the Simojoki river basin in northern Finland and at Dalelva Brook in northern Norway. The new modules did not imply any significant changes in simulated NO3- concentration levels in the streams but improved the timing of simulated higher concentrations. The inclusion of a modified temperature response function and an empirical snow-correction factor improved the flexibility and applicability of the model for climate effect studies.

Heat tolerance in groundnut

Tolerance to high soil and air temperature during the reproductive phase is an important component of adaptation to and and semi-arid cropping environments in groundnut. Between 10 and 22 genotypes were screened for tolerance to high air and soil temperature in controlled environments. To assess tolerance to high soil temperature, 10 genotypes were grown from start of podding to harvest at ambient (28 degrees) and high (38 degreesC) soil temperatures, and crop growth rate (CGR), pod growth rate (PGR) and partitioning (ratio PGR:CGR) measured. To assess tolerance to high air temperature during two key stages-microsporogenesis (3-6 days before flowering, DBF) and flowering, fruit-set was measured in two experiments. In the first experiment, 12 genotypes were exposed to short (3-6 days) episodes of high (38 degreesC) day air temperature at 6 DBF and at flowering. In the second experiment, 22 genotypes were exposed to 40 degreesC day air temperature for I day at 6 DBF, 3 DBF or at flowering. Cellular membrane thermostability (relative injury, RI) was also measured in these 22 genotypes. There was considerable variation among genotypes in response to high temperature, whether assessed by growth rates, fruit-set or RI. Pod weight at high soil temperature was associated with variation in CGR rather than partitioning. Flowering was more sensitive to high air temperature than microsporogenesis. Genotypes tolerant to high air temperature at microsporogenesis were not necessarily tolerant at flowering, and nor was tolerance correlated with RI. Six genotypes (796, 55-437, ICG 1236, ICGV 86021, lCGV 87281 and ICGV 92121) were identified as heat tolerant based on their performance in all tests. These experiments have shown that groundnut genotypes can be easily screened for reproductive tolerance to high air and soil temperature and that several sources of heat tolerance are available in groundnut germplasm. (C) 2003 Elsevier Science B.V. All rights reserved.

Drought, pod yield, pre-harvest Aspergillus infection and aflatoxin contamination on peanut in Niger

Soil moisture and soil temperature affect pre-harvest infection with Aspergillus flavus and production of aflatoxin. The objectives of our field research in Niger, West Africa, were to: (i) examine the effects of sowing date and irrigation treatments on pod yield, infection with A. flavus and aflatoxin concentration; and (ii) to quantify relations between infection, aflatoxin concentration and soil moisture stress. Seed of an aflatoxin susceptible peanut cv. JL24 was sown at two to four different sowing dates under four irrigation treatments (rainfed and irrigation at 7, 14 and 21 days intervals) between 1991 and 1994, giving 40 different 'environments'. Average air and soil temperatures of 28-34 degrees C were favourable for aflatoxin contamination. CROPGRO-peanut model was used to simulate the occurrence of moisture stress. The model was able to simulate yields of peanut well over the 40 environments (r(2) = 0.67). In general, early sowing produced greater pod yields, as well as less infection and lower aflatoxin concentration. There were negative linear relations between infection (r(2) = 0.62) and the average simulated fraction of extractable soil water (FESW) between flowering and harvest, and between aflatoxin concentration (r(2) = 0.54) and FESW in the last 25 days of pod-filling. This field study confirms that infection and aflatoxin concentration in peanut can be related to the occurrence of soil moisture stress during pod-filling when soil temperatures are near optimal for A. flavus. These relations could form the basis of a decision-support system to predict the risk of aflatoxin contamination in peanuts in similar environments. (c) 2005 Elsevier B.V. All rights reserved.

The effect of wheat straw residue on the emergence and early growth of sugar beet (Beta vulgaris) and oilseed rape (Brassica napus)

The management of straw residue can be a concern in non-inversion tillage systems where straw tends to be incorporated at shallow depths or left on the soil surface. This can lead to poor crop establishment because straw residue can impede or hinder crop emergence and growth. Small container-based experiments were undertaken using varying amounts of wheat straw residue either incorporated or placed oil the soil surface. The effects on (lays to seedling emergence, percentage emergence, seedling dry-weight and soil temperature using sugar beet and oilseed rape were investigated because these crops often follow wheat in a cropping sequence. The position of the straw residue was found to be the primary factor in reducing crop emergence and growth. Increasing the amount of straw residue (from 3.3 t ha(-1) to 6.7 t ha(-1)) did not show any consistent trends in reducing crop emergence or growth. However, in some instances, results indicated that an interaction between the position and the amount of straw residue Occurred particularly when the straw and seed was placed on the soil surface. Straw placed on the soil surface significantly reduced mean day-time soil temperature by approximately 2.5 degrees C compared to no residue. When the seed and straw was placed on the soil Surface a lack of seed-to-soil contact caused a reduction in emergence by approximately 30% because of the restriction in available moisture that limited the ability for seed imbibition. This trend was reversed when the seed was placed in the soil, but with straw residue still on the soil surface, because the surface straw was likely to reduce moisture evaporation and improved seed-to-soil contact that led to rapid emergence. In general, when straw was mixed in or placed on the soil surface along with the seed, sugar beet and oilseed rape emergence and early growth biomass was significantly restricted by approximately 50% compared to no residue. The consequences of placing seed with or near to straw residue have been shown to cause a restriction in crop establishment. In both oilseed tape and sugar beet, this could lead to a reduction in final crop densities, poor, uneven growth and potentially lower yields that could lower financial margins. Therefore, if farmers are planning to use non-inversion tillage methods for crop establishment, the management and removal of straw residue from near or above the seed is considered important for successful crop establishment. (C) 2008 Elsevier B.V. All rights reserved.

Influence of drought-induced acidification on the mobility of dissolved organic carbon in peat soils

A strong relationship between dissolved organic carbon (DOC) and sulphate (SO42−) dynamics under drought conditions has been revealed from analysis of a 10-year time series (1993–2002). Soil solution from a blanket peat at 10 cm depth and stream water were collected at biweekly and weekly intervals, respectively, by the Environmental Change Network at Moor House-Upper Teesdale National Nature Reserve in the North Pennine uplands of Britain. DOC concentrations in soil solution and stream water were closely coupled, displaying a strong seasonal cycle with lowest concentrations in early spring and highest in late summer/early autumn. Soil solution DOC correlated strongly with seasonal variations in soil temperature at the same depth 4-weeks prior to sampling. Deviation from this relationship was seen, however, in years with significant water table drawdown (>−25 cm), such that DOC concentrations were up to 60% lower than expected. Periods of drought also resulted in the release of SO42−, because of the oxidation of inorganic/organic sulphur stored in the peat, which was accompanied by a decrease in pH and increase in ionic strength. As both pH and ionic strength are known to control the solubility of DOC, inclusion of a function to account for DOC suppression because of drought-induced acidification accounted for more of the variability of DOC in soil solution (R2=0.81) than temperature alone (R2=0.58). This statistical model of peat soil solution DOC at 10 cm depth was extended to reproduce 74% of the variation in stream DOC over this period. Analysis of annual budgets showed that the soil was the main source of SO42− during droughts, while atmospheric deposition was the main source in other years. Mass balance calculations also showed that most of the DOC originated from the peat. The DOC flux was also lower in the drought years of 1994 and 1995, reflecting low DOC concentrations in soil and stream water. The analysis presented in this paper suggests that lower concentrations of DOC in both soil and stream waters during drought years can be explained in terms of drought-induced acidification. As future climate change scenarios suggest an increase in the magnitude and frequency of drought events, these results imply potential for a related increase in DOC suppression by episodic acidification.

Five years of carbon dioxide fluxes measurements in a highly vegetated suburban area

Suburban areas continue to grow rapidly and are potentially an important land-use category for anthropogenic carbon-dioxide (CO2) emissions. Here eddy covariance techniques are used to obtain ecosystem-scale measurements of CO2 fluxes (FC) from a suburban area of Baltimore, Maryland, USA (2002–2006). These are among the first multi-year measurements of FC in a suburban area. The study area is characterized by low population density (1500 inhabitants km−2) and abundant vegetation (67.4% vegetation land-cover). FC is correlated with photosynthetic active radiation (PAR), soil temperature, and wind direction. Missing hourly FC is gap-filled using empirical relations between FC, PAR, and soil temperature. Diurnal patterns show net CO2 emissions to the atmosphere during winter and net CO2 uptake by the surface during summer daytime hours (summer daily total is −1.25 g C m−2 d−1). Despite the large amount of vegetation the suburban area is a net CO2 source of 361 g C m−2 y−1 on average.

Evidence of increased net ecosystem productivity associated with a longer vegetated season in a deciduous forest in south-central Indiana, USA

Observations of net ecosystem exchange (NEE) of carbon and its biophysical drivers have been collected at the AmeriFlux site in the Morgan-Monroe State Forest (MMSF) in Indiana, USA since 1998. Thus, this is one of the few deciduous forest sites in the world, where a decadal analysis on net ecosystem productivity (NEP) trends is possible. Despite the large interannual variability in NEP, the observations show a significant increase in forest productivity over the past 10 years (by an annual increment of about 10 g C m−2 yr−1). There is evidence that this trend can be explained by longer vegetative seasons, caused by extension of the vegetative activity in the fall. Both phenological and flux observations indicate that the vegetative season extended later in the fall with an increase in length of about 3 days yr−1 for the past 10 years. However, these changes are responsible for only 50% of the total annual gain in forest productivity in the past decade. A negative trend in air and soil temperature during the winter months may explain an equivalent increase in NEP through a decrease in ecosystem respiration.

An energy budget agent-based model of earthworm populations and its application to study the effects of pesticides

Earthworms are important organisms in soil communities and so are used as model organisms in environmental risk assessments of chemicals. However current risk assessments of soil invertebrates are based on short-term laboratory studies, of limited ecological relevance, supplemented if necessary by site-specific field trials, which sometimes are challenging to apply across the whole agricultural landscape. Here, we investigate whether population responses to environmental stressors and pesticide exposure can be accurately predicted by combining energy budget and agent-based models (ABMs), based on knowledge of how individuals respond to their local circumstances. A simple energy budget model was implemented within each earthworm Eisenia fetida in the ABM, based on a priori parameter estimates. From broadly accepted physiological principles, simple algorithms specify how energy acquisition and expenditure drive life cycle processes. Each individual allocates energy between maintenance, growth and/or reproduction under varying conditions of food density, soil temperature and soil moisture. When simulating published experiments, good model fits were obtained to experimental data on individual growth, reproduction and starvation. Using the energy budget model as a platform we developed methods to identify which of the physiological parameters in the energy budget model (rates of ingestion, maintenance, growth or reproduction) are primarily affected by pesticide applications, producing four hypotheses about how toxicity acts. We tested these hypotheses by comparing model outputs with published toxicity data on the effects of copper oxychloride and chlorpyrifos on E. fetida. Both growth and reproduction were directly affected in experiments in which sufficient food was provided, whilst maintenance was targeted under food limitation. Although we only incorporate toxic effects at the individual level we show how ABMs can readily extrapolate to larger scales by providing good model fits to field population data. The ability of the presented model to fit the available field and laboratory data for E. fetida demonstrates the promise of the agent-based approach in ecology, by showing how biological knowledge can be used to make ecological inferences. Further work is required to extend the approach to populations of more ecologically relevant species studied at the field scale. Such a model could help extrapolate from laboratory to field conditions and from one set of field conditions to another or from species to species.

ERA-Interim/Land: a global land surface reanalysis data set

ERA-Interim/Land is a global land surface reanalysis data set covering the period 1979–2010. It describes the evolution of soil moisture, soil temperature and snowpack. ERA-Interim/Land is the result of a single 32-year simulation with the latest ECMWF (European Centre for Medium-Range Weather Forecasts) land surface model driven by meteorological forcing from the ERA-Interim atmospheric reanalysis and precipitation adjustments based on monthly GPCP v2.1 (Global Precipitation Climatology Project). The horizontal resolution is about 80 km and the time frequency is 3-hourly. ERA-Interim/Land includes a number of parameterization improvements in the land surface scheme with respect to the original ERA-Interim data set, which makes it more suitable for climate studies involving land water resources. The quality of ERA-Interim/Land is assessed by comparing with ground-based and remote sensing observations. In particular, estimates of soil moisture, snow depth, surface albedo, turbulent latent and sensible fluxes, and river discharges are verified against a large number of site measurements. ERA-Interim/Land provides a global integrated and coherent estimate of soil moisture and snow water equivalent, which can also be used for the initialization of numerical weather prediction and climate models.